The higher order C n dispersion coefficients for 
the alkali atoms. 
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The van der Waals coefficients, from Cn through to Cie resulting from 2nd, 3rd and 4th order 
perturbation theory are estimated for the alkali (Li, Na, K and Rb) atoms. The dispersion co- 
efficients are also computed for all possible combinations of the alkali atoms and hydrogen. The 
parameters are determined from sum-rules after diagonalizing the fixed core Hamiltonian in a large 
basis. Comparisons of the radial dependence of the C„/r n potentials give guidance as to the radial 
regions in which the various higher-order terms can be neglected, ft is seen that including terms up 
to Cio/r 10 results in a dispersion interaction that is accurate to better than 1% whenever the inter- 
nuclear spacing is larger than 20 ao. This level of accuracy is mainly achieved due to the fortuitous 
cancellation between the repulsive (Cn, C13, C15) and attractive (C12, C14, C\e) dispersion forces. 

PACS numbers: 34.20.Cf, 31.25.Jf, 31.15.Pf, 32.10.Dk 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sation (BEC) for the alkali atoms, Li, Na, Rb and atomic 
hydrogen [lj, |2j has resulted in an upsurge of interest in 
the area of cold atom physics. One consequence of this is 
the increased importance in determining the interaction 
potentials between alkali and alkaline-earth atoms. For 
example, the stability and structure of BECs depends on 
the sign (and magnitude) of the scattering length, and 
the scattering length sensitively on the fine details of the 
interaction potential 0, 0] . 

One can try and determine the interaction potential by 
either explicit calculation 0, 0, Q , or by analysis of high 
precision spectroscopic experiments of the dimer [H ]j| 
LfJ. These approaches are most suited to determining the 
potential at small to intermediate values of internuclear 
distance (e.g. r < 25 ao). At large distances, one relies 
on the approach pioneered by Dalgarno and collaborators 
[ill Il2l Il3|| , namely using oscillator strength sum rules 
to determine the so-called dispersion coefficients. 

The long-range interaction between two spherically 
symmetric atoms can be written in the general form 

V(r - 00) = -V 9 (r)-V s {r)-V 10 (r)-Vu(r)-V 12 (r)-. . . , 
where the dispersion potential, V„ (r) of order n is written 



(2) 



The C n parameters are the van der Waals dispersion co- 
efficients. The even (n = 6, 8, . . .) dispersion coefficients 
are calculated using oscillator-strength sum rules derived 
from 2nd-order perturbation theory and provide an at- 
tractive interaction. The odd (n = 11, 13, . . .) terms 
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come from 3rd-order perturbation theory, and are re- 
pulsive |14| . Contributions from 4th-order perturbation 
theory start at n = 12 [l4l HEj. The calculations were 
stopped at n = 16 since the contributions from 5th and 
6th-order perturbation theory start at n — 17 and n = 18 
respectively. 

It has been customary in previous systematic studies of 
dispersion coefficients to restrict the calculations to just 
the Cq, C$ and C10 terms 16]. With the recent inter- 
est in determining the scatteri ng le ngths of A-A systems 
to high precision 0, 0, 0, 0, |2(j , it is now worth- 
while to evaluate the higher n dispersion coefficients with 
a view to determining whether these make a significant 
contribution to the inter-atomic potential for those radial 
separations at which the long-range form of the potential 
is first applied. 

The recent spectrum analyses have typically an inter- 
nuclear distance of about 20 ao as the boundary point 
used to join the explicitly determined potential with the 
asymptotic form given in terms of the dispersion forces 0, 
Taking the Rb-Rb potential as a specific example, 
explicit calculations have shown that Cq — 4.635 x 10 3 
au, C* 8 = 5.701 x 10 5 au, C w = 7.916 x 10 7 au, and 
C12 = 1.427 x 10 10 au 21). At a radius of 20 a the 
ratios of the contributions from the various terms are 1 : 
0.308 : 0.1068 : 0.0481. The contribution of the highest 
multipole potential is about 5% of the dipole-dipole Cq 
term. 

In this article, the Cn to C\q dispersion coefficients 
and related parameters are computed using a semi- 
empirical approach for the alkali atoms. The method 
adopted utilizes oscillator strengths calculated within a 
semi-empirical framework |2ll . Comparisons with previ- 
ous high accuracy ab-initio calculations [2ll I22I I23L I24I 
l25| of the Cq , C$ and C10 dispersion coefficients suggests 
the method is able to generate coefficients that are accu- 
rate to about 1-2%. 
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II. CALCULATION OF THE DISPERSION 
PARAMETERS 

All the dispersion coefficients computed in this paper 
were computed by first diagonalizing the semi-empirical 
Hamiltonian in a large mixed Laguerre type orbital or 
Slater type orbital basis set. Next various sum rules in- 
volving oscillator strengths or reduced matrix elements 
were summed over the set of physical and pseudo states. 

The respective formulae for the various coefficients are 
now given. The notation C n is used to denote the total 
dispersion coefficient for a given n. So for n > 12, C n 
is the sum of C {2) and Ci 4) where c£ 2) arises from 2nd- 
order perturbation theory and arises from 4th-order 
perturbation theory. 



A. The 2nd order contributions to Cn and C14 



The sum rules involve contributions from both core 
and valence excitations. The valence contributions were 
evaluated by diagonalizing the model Hamiltonian in a 
very large basis. Determination of the /-value distribu- 
tion for the core was handled by using the properties of 
/-value sum rules and an approach that gives a reason- 
able estimate of the /-value distribution with a minimum 
of computation 21]. We use the sum-rule for the polar- 
izability, eq. (@J and 



£N(: 



,2l-2\ 



> = £/f=^)(o) 



(6) 



|2lll31| to estimate an /^-value distribution function of 
reasonable accuracy. This expression reduces to the well 
known Thomas- Reiche-Kuhn sum rule for 1=1, viz 



N 



£/f } = s (1) (o) 



(7) 



While the higher n dispersion coefficients for the 
H-H system have been determined to high precision 
|26l l27i l28j , only a few calculations of C\i and C14 have 
been done for the alkali atoms 0, |2|j . The calculations 
of Ovsiannikov, Guilyarovski and Lopatko (OGL) [T4II2II 
used an approximate expression for the Greens function 
and can only be expected to be accurate at the 20% level. 
Explicit comparisons with the C% and Cm values of Patil 
and Tang (PT) [2J| in references [21], |25j reveal discrep- 
ancies of order 10% for Na2, K 2 and Rb 2 . Recursion rules 
exist for estimating C n at high n from data tabulations 
of C n for lower values of n |27fl . 

The polarization and dispersion parameters can be 
computed from their respective oscillator strength sum- 
rules which are well known. The oscillator strength, 
from the ground state (with orbital and spin angular mo- 
mentum equal zero) to the ith excited state is defined as 



In these expressions N is the total number of electrons 



{i) = 2|(^q || r*C\v) || ^)| 2 e 0l 
Ioi (2£+l) 



(3) 



In this expression C e is the spherical tensor of rank £, the 
Wigner-Eckart theorem is defined according to jsjj, and 
eoi is the excitation energy of the transition. The sum 
rule for the adiabatic multipole polarizability, a''' is 



4 6 0i 



(4) 



The 2nd-order contributions to C12 and C14 are defined 
, (2) _ ^ (2n-2)! 



f {U) f lj) 



2 " 2^-! ^ eoieo .( eoi+eoi ) 



(5) 



with ii + £j + 1 = n. The sum rules are a generalized sum 
which implicitly includes a sum over excitations to bound 
states and an integration taking into account excitations 
to continuum states. 



and 



is an expectation value of the ground state 



wave function. 

The core /(^ -value distribution was determined by 
assuming each closed sub-shell made a contribution of 
N ir f- 2 to the sum-rules, where Ni the number of elec- 
trons in the sub-shell and r 2l ~ 2 is the ( r 2t ~ 2 ) Hartree- 
Fock (HF) expectation value. Next, the excitation energy 
for each sub-shell is set to the HF single particle, £j plus 
an energy shift. The energy shift, was set by using 

(£) 

the core multipole polarizability, acore and the relation 



a. 



E 



21-2 



(e; + AW)2 



(8) 



A full description of the details and core polarizabilitics 
used to fix A^) for I = 1,2 and 3 has been published [2l| . 
The contribution from the core was omitted for £ = 4 and 
5 since there are no reliable estimates of these core po- 

(2£— 2) 

larizabilities and the £r\ weighting factors lead to 
the core becoming less important for larger £. The core 
contribution to was less than 0.2% for Rb and most 
of the core contribution to comes from the /W dis- 
tribution. 



B. The 3rd order and 4th order potentials 



The dispersion coefficients, Cji and C p, a rise from 
3rd-order perturbation theory [3 Q HI M, HI 
The standard expressions are expressed in terms of sums 
of products of reduced matrix elements. The expressions 
derived by OGL 0111 are used. Given the complexity 
of the expressions, it is not surprising there have been 
relatively few calculations of these coefficients for any 
atoms. Accurate values for Cn through to Qji for the 
H-H system have recently been published [2H . |28| . 

For the alkali atoms, the only estimates of Cn have 
been obtained from the relatively small calculations of 
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OGL 14] and Patil and Tang [29I There have been 
some Cii and C\2 estimates for Cs [33, but that atom is 
deemed too heavy to accurately describe with the present 
non-relativistic method and is not discussed any further 
here. 

The dispersion coefficients, C\% and C14, both have 
contributions that arise from 4th-order perturbation the- 
ory. However, there have been very few calculations of 
the 4th order term for any systems. Bukta and Meath 
(llj made an explicit calculation f° r the H-H dimer. 
Besides deriving a general expression OGL |14l |26J also 

made some estimates of C\ 2 for combinations of hydro- 
gen and the alkali atoms. Finally, the present authors 
made a comprehensive calculation of c4 up to n = 32 
for the H-H interaction [2(| using the OGL formalism. 

Although the core has not been taken explicitly into 
account in the evaluation of the 3rd and higher order dis- 
persion parameters their impact is expected to be smaller 
than 2nd order due since the 3rd and 4th order sum-rules 
have an energy denominator involving the square or cube 
of the excitation energy. 



TABLE I: The higher multipole polarizabilities for the lighter 
alkali atoms. All values are in atomic units. 



Method 


1CT 6 


Q< 4 ' 


lO" 9 


10-9 (6) 




present 


PT [291 


Present 


Present 


Li 


1.999 


1.947 


0.1551 


17.04 


Na 


2.973 


2.828 


0.2450 


28.53 


K 


11.82 


10.69 


1.204 


172.4 


Rb 


16.57 


14.49 


1.762 


259.6 



III. RESULTS OF THE CALCULATION 

A. Sensitivity of calculations to the ground state 
wave function 

One aspect of the calculations that warrants particular 
mention was the sensitivity of the higher-^ polarizabili- 
ties (listed in Tabled) and thus C$ , c[f and C<$ to the 
representation of the ground state wave function. While 
there are many advantages to representing the wave func- 
tion with a linear combination of convenient basis func- 
tions, there are some negative features. One of those neg- 
ative features relates to the behavior at large distances 
from the nucleus. Unlike a grid based calculation, the 
correct asymptotics are not imposed and so the large-r 
part of the wave function, which has a weak influence on 
the total binding energy can be inaccurate. 

This is best illustrated by a specific calculation. Our 
initial calculations for the and polarizibilities Na 
used a 3s wave function written as a linear combination 
of 12 Slater type orbitals (STOs). This wave function had 
a binding energy against ionization of 0.1888532 hartree. 



The resulting polarizabilities were offi = 3.46 x 10 6 au 
and a (5) = 6.42 x 10 s au. When the STO basis was 
replaced by a large Laguerre type orbital (LTO) basis 
(this was necessary for the evaluation of Cn) and the 
energy driven to convergence the resulting binding en- 
ergy was 0.1888549 hartree. However there were dra- 
matic changes in the polarizabilities, with the new val- 
ues being = 2.97 x 10 6 au and = 2.45 x 10 8 au 
(refer to Table [IJ. The polarizability decreased by 
a factor of about 2.5 when the binding energy changed 

by 1.7 x 10~ 6 hartree! The dispersion parameters, 

and (7^2 and were also sensitive to the representation 

(2) 

of the ground state wave function. The initial C{ 4 of 
4.04 x 10 11 a.u. was decreased to 2.602 x 10 11 a.u. when 
the ground state basis was made exhaustively large. The 
impact of the basis set on the dispersion coefficients was 
was not so extreme as for the polarizability but was still 
substantial. 

The /( 2 ) and /( 3 ) distributions used in earlier calcu- 
lations of the Na and K polarizabilities and dispersion 
coefficients 0, |^ were affected to a smaller extent by 
this problem. The octupole polarizabilities were 3-4% 
too large for these atoms, while the effect upon C10 was 
to make them about 1-2% too large. The impact upon 
o^ 2 -* and Cg was an order of magnitude smaller but it 
was still discernible [2(| . Where values of Cs and Cjn are 
required in the present work, the revised values |35j are 
used. 



B. The homo- nuclear Alkali- Alkali systems 

The results of the present calculations for the Li, Na, 
K and Rb homo-nuclear alkali atom pairs are listed in 
Tables [I] and [HJ and compared with the calculations of 
Ovsiannikov et al and Patil and Tang [2{j . 

(2) 

All of the present values of C\ 2 are larger than the PT 
and OGL estimates by amounts ranging from 10-50%. 
This is not a concern since the OGL and PT estimates 
of Cq, Cs and Cm are also smaller than the latest data 
for these parameters 0, 0, HH ■ The present estimates 
of Cn and C13 are also 10-50% larger than the PT and 
OGL estimates. 

In TableHH ci 2) and ci 4) are given as separate entries. 
The total dispersion parameter C n is given in Table IIVI 
which tabulates the dispersion parameters for all possible 

atom-atom combinations. The 4th-order terms, C^ are 

(2) 

about 4-5% the size of C 12 . The 4th-order correction 
to C12 is about the same size as the correction due to 
the core. The size of the 4th-order correction to C14 is 
larger, with about 25% of the final value of C14 coming 
from €'[2' '• F° r C\§, the 4th-order terms are almost as 
large as the 2nd-order terms, with C^ being 80% the 
size of C$ for K2 and Rb2- 

The contribution of the core to , C$ and c[^ is 
relatively small. The contribution is largest for the Rb2 
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(2) 

TABLE II: The dispersion coefficients for homo-nuclear combinations of the lighter alkali atoms. The C„ terms have core 
contributions. All values are in atomic units. 

Method 1Q- 6 Cn 1(T 9 Cff 1(T 9 Cff 1Q- 9 C 13 HT 12 C™ HT 12 Cff HT 12 C 15 HT 15 C™ HT 15 C[f 

Li 

Present -40.44 0.9015 0.0401 -11.05 0.1455 0.04251 -2.873 0.02997 0.01908 

PT [22] -37.36 0.8648 -10.11 

OGL [14] -36.0 0.648 0.0433 -8.85 

Na 

Present -61.01 1.500 0.06699 -18.90 0.2602 0.0720 -5.368 0.05736 0.03637 

PT [22] -53.32 1.375 -16.46 

OGL [14] -39.5 0.917 0.0396 -11.3 

K 

Present -364.9 9.102 0.4640 -147.2 2.000 0.6807 -53.02 0.5544 0.4480 

PT [22] -312.3 7.749 -122.8 

OGL [14] -250 5.430 0.361 -89.3 

Rb 

Present -536.2 14.26 0.6986 -236.5 3.314 1.098 -91.20 0.9656 0.7897 

PT [29] -464.4 11.56 -195.1 

OGL [14] -376 8.370 0.556 -145.0 



(2) 

dimer, but even here the effect is 4.2% for C\t> , 3.2% for 
d 2 ? and 2.5% for C$ . 



C. Critical radius for dispersion formula 

The LeRoy radius is often used as an estimate of the 
critical radius beyond which the interaction can be de- 
scribed by the use of a dispersion formula [3(J |37j . It is 
defined for two atoms, A and B as 

Rlr^2(^U + ^)b) , (9) 

where (r 2 )^ is evaluated for the ground state. The LeRoy 
radius for the homo-nuclear dimers are given in Table ITTT1 
We introduce some parameters so that the range of 
validity of the dispersion formula can be discussed in a 
quantitative manner. First, the partial sum of the dis- 
persion energy up to the nth term is defined as 

ra—n 

W n (r) = V m(r) • (10) 

m— 6 

One would then consider the relative size of W n (r) to 
Woo(r) as a measure of the accuracy of the truncated 
dispersion potential to the exact potential. There are of 
course problems associated with the evaluation of Woo, 
and so we identify Woo(r) with Wxe(r). It is natural to 
stop the analysis at n = 16 since the contributions from 
5th and 6th-order perturbation theory start at n = 17 
and n = 18 respectively. 

FigureGJshows the ratio |Wi.o(r) — W r i 6 (r)|/|W / i 6 (r)| as 
a function of r for the Li, Na, K and Rb dimers. This ra- 
tio is seen to be smaller than 1% for all values of r greater 



than 20 arj. The curves for Li2 and Na2 are very simi- 
lar, as are the curves for K2 and Rb2. The existence of 
two sets of two very similar curves is probably related to 
the fact that the lowest lying d-excitations for Li and Na 
involve a change in principal quantum number whereas 
those for K and Rb do not. 

There is a sign change in (Wio(r) — Wie(r)) near 20 ao 
for all four dimers. The magnitude of Wio(r) is generally 
smaller than W\§(r) for small separations, but for large 
separations W\q (r) is generally larger in magnitude than 
W\e(r). This is caused by the repulsive Vn(r), Vis(r) 
and Vi (r) interactions. 



TABLE III: Various radial distances (in ao) related to the 
accuracy of different order expansions of the dispersion pa- 
rameters for the homo-nuclear alkali-metal atom dimer. 



Atom 


<r 2 > 


Rlr 


Re 


Rg 


Rio 


Li 


17.47 


16.72 


77.53 


25.92 


16.30 


Na 


19.51 


17.67 


86.24 


28.07 


17.24 


K 


27.97 


21.16 


103.8 


32.61 


18.82 


Rb 


30.76 


22.18 


111.0 


34.44 


19.54 



A useful way to parameterize this information is to 
define R n such that it gives the smallest radius for which 
the partial dispersion energy W n (r) is accurate to 1%. In 
effect, R n is the largest r solution of the equation 

\W n (r) - Wx 6 (r)\ = 0.010 x W w {r) n < 16 (11) 

Table IIIII gives R n for the hydrogen and alkali homo- 
nuclear dimers. The values of R n get larger as the atom 
gets heavier. A dispersion potential only involving C§ 
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would not be accurate to 1% until the separation dis- 
tance increased to more than 100 ao in the case of Rb 2 . 
This distance shrinks dramatically with the inclusion of 
the Vs and Vio potentials. The size of Rs indicates that 
a dispersion interaction with only Vq and V% is not good 
enough to describe these alkali dimers. The R\q parame- 
ter is smaller than the LeRoy radius for all systems, and 
so the use of W\o{r) will be accurate to better than 1% as 
long as the internuclear separation is greater than Rlr- 

FiguresEHHlHandElshowthe ratio |K(r-)|/|Wi 6 (r)| as 
a function of r for all the alkali dimers. One of the most 
noticeable features of these curves is the existence of a 
nexus point for the Vn, V13, V14, V15 and Viq potentials 
between 16-21 a for all alkali dimers. This suggests that 
a recursion relation of the type C n + m — A m C n exists 
between the respective dispersion coefficients. 

The biggest correction to Wi${r) for sufficiently large 
r should be V\\{r). The C\2-C\\ ratios for Li 2 , Na 2 , 
K 2 and Rb 2 from Table llVl are 23.3, 25.7, 26.2 and 27.9 
respectively and, for, r less than these values it is actually 
Vi2(r) which is the largest correction to W\o(r) (note; at 
sufficiently small r, V13, V14, . . . will also be larger than 
Vxi). 

The extent to which mutual cancellations act to mini- 
mize the error can be gauged by adding the magnitudes 
of V n (r) at r = 20. For Na, the J^HLn l^«( r )l i s equal 
to 0.047 x Wie(r) at r = 20a . But, the partial sum 
|W 10 (r) - W 16 (r)\ is equal to 0.0012 x W 16 (r) at this 
radius. The alternating signs reduce the impact of the 
higher order terms by a factor of about 40. For Rb, 
the E^Lii \V n (r)\ sum gives 0.168 x W 16 {r) at r = 20 
clq. The magnitudes of |Ki(r)|/|Wio(r-)| are greater than 
0.025 for all n between 11 and 16. The partial sum 
\W w (r) - W w (r)\ is equal to 0.0071 x W 16 (r) at this 
radius. In this case, the alternating signs have reduced 
the impact of the higher order terms by a factor of more 
than 20. 

Figures El and [7| better illustrate the extent to which 
mutual cancellations result in W\o being a very good ap- 
proximation to the total dispersion potential for Na and 
Rb. The combination (Vn(r) + Vi 2 (r)) is much smaller 
than either Vn or V12, and similarly (Vi3(r) + Vu(r)) 
is smaller than V13 or V14 while (Vis(r) + Vi§(r)) is 
smaller than V\§(r) or V\e{r). The combination of 
V 2n -i(r) + V 2n (r) falls to 1% of W w (r) about 5-10 a 
closer to the origin than either V 2 „_i(r) or V 2rl (r). 



1. The Rubidium dimer 

A recent analysis of three complementary experiments 
by van Kempen et al for the Rb 2 dimer resulted in 
experimental estimates of the dispersion parameters Cq, 
Cg, C10 and C\\ and gave estimates of the scattering 
length to a precision of about 1% for Rb 87 . 

The present work suggests that their estimates of Cq, 
C$ and C10 may be well founded and our previous cal- 
culations agreed with the experimental values to an ac- 




r (units of a ) 

FIG. 1: Theratioof \Wio(r)-Wi 6 (r)\/\Wi 6 (r)\ asafunction 
of r (units of ao) for all the possible homo-nuclear pairs. 

curacy of 2% 21] . However, it is probable that their 
attempt to determine Cn was overly ambitious. The 
present calculation gave Cn = —5.362 x 10 9 au which 
lies outside the van Kempen estimate of Cn = (—8.6 ± 
0.17) x 10 9 au. We believe this difference is a consequence 
of assumptions made by van Kempen et al when they 
performed the fit to extract the dispersion coefficient. 

First, their least squares fit relied on a value of Ci 2 of 
1.19 x 10 10 au computed by Patil and Tang 38]. This 
estimate of Ci 2 is about 25% smaller than the present 
2nd plus 4th order Ci 2 of 1.496 x 10 10 au. Second, their 
estimate is derived from the energies of bound states that 
have an outer turning radius of about 20 ao. At this ra- 
dius, V\2 is actually larger than Vn. Furthermore, the 
higher order terms, V13, V14, F15 and Vie were all slightly 
larger than Vn at a radius of 20 ao! However, as has 
been noted, the alternating signs of the successive terms 
results in a considerable degree of cancellation. Never- 
theless, one can conclude that the van Kempen estimate 
of Cn is sensitive to the accuracy and presence of the 
high n terms in the dispersion interaction. 

D. The hetero- nuclear systems including H 

The n > 10 dispersion parameters for all possible com- 
binations of H, Li, Na, K and Rb are given in Table HVI 
The radial matrix elements for hydrogen were those of 
the N = 15 calculation used in an earlier calculation of 
the H-H dispersion parameters [2(|. TheCSi and C\2 pa- 
rameters reported by OGL 0] and PT |29j are not listed 
since the present calculations are more sophisticated. 

IV. CONCLUSIONS 

The complete set of dispersion parameters up to n = 16 
has been computed for all combinations of hydrogen and 
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TABLE IV: The dispersion coefficients, Cn — > Cie for all the possible interacting pairs formed by hydrogen and the alkali 
atoms. The C12 , C14 and Cie coefficients have contributions from both 2nd- and 4th-order perturbation theory. All values are 
in atomic units. 



System 


C11 


C12 


C*13 


C14 


C*15 




H-H [26] 


-3474.9 


1.2273 x 10 5 


-3.2699 x 10 5 


6.3617 x 10 6 


-2.8396 x 10 7 


4.4121 x 10 s 




10 -6 Cu 


10 -6 C12 


10 -9 C13 


10 -9 C14 


10 -9 C15 


10 -12 Cie 


H-Li 


-0.2251 


21.17 


-0 04752 


2.735 


-9.438 


0.4593 


H-Na 


-0.2760 


29.43 


-0 06269 


4.030 


-13.25 


0.7160 


H-K 


-0 5636 


92.69 


-0.1635 


15.68 


-43.10 


3.409 


H-Rb 


-0.6619 


12.34 


-0.2064 


21.86 


-57.57 


4.936 




10 -6 C11 


10 -9 Cl9 


10 -9 C^% 

-LVJ w 13 


10~ 12 C14 


10 -12 CV 


10 -13 Cie 


Li-Li 


-40.44 


0.9417 


-11.05 


0.1880 


-2.873 


4.906 


Li-Na 


-49.63 


1.221 


-14.48 


0.2512 


-3.937 


6.819 


Li-K 


-119.4 


3.216 


-40.67 


0.7685 


-12.68 


24.08 


Li-Rb 


-143.9 


4.134 


-51.84 


1.018 


-16.84 


32.90 


Na-Na 


-61.01 


1.567 


-18.90 


0.3323 


-5.368 


9.373 


Na-K 


-146.5 


4.015 


-52.60 


0.9834 


-17.04 


31.97 


Na-Rb 


-176.5 


5.125 


-66.85 


1.292 


-22.54 


43.32 


K-K 


-364.9 


9.567 


-147.2 


2.681 


-53.02 


100.24 


K-Rb 


-442.1 


11.99 


-186.7 


3.449 


-69.63 


133.0 


Rb-Rb 


-536.2 


14.96 


-236.5 


4.412 


-91.20 


175.5 



10 u F 
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FIG. 2: The ratio of |Ki(f)|/|W / i6(r)| as a function of r (units 
of ao) for the lithium dimer. The decay of the ratio at large 
r is fastest for large n. 
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FIG. 3: The ratio of |K(r)|/|Wi 6 (r)| as a function of r (units 
of ao) for the sodium dimer. The decay of the ratio at large 
r is fastest for large n. 



the alkali atoms up to rubidium. The relative importance 
of the dispersion potentials V n (r) increases as the atoms 
get heavier. It was found that the dispersion energy given 
by the first 3 terms of eq. is accurate to 1% whenever 
R > 20 ao- This degree of accuracy at relatively small 
internuclear separations comes from a fortuitous cancel- 
lation between the terms with n > 10 in the dispersion 
energy. The 3rd-order potentials, Vn, V43 and V45 are 



repulsive, while the even terms, Via, Vu an d Vl6 are at- 
tractive. The 4th-order contribution to Cie was almost 
as large as the 2nd-order contribution to Ciq . 

Whether terms in the dispersion interaction with n > 
10 are important in the description of alkali dimers is 
essentially a question about whether the dispersion in- 
teraction has to be known to a precision of 1.0% or 0.1%. 
It is clear that additional terms going beyond C10 should 
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FIG. 4: The ratio of |Ki(r)|/|VKi6(V)| as a function of r (units 
of ao) for the potassium dimer. The decay of the ratio at large 
r is fastest for large n. 
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FIG. 5: The ratio of |Ki(f)|/|W / i6(r)| as a function of r (units 
of ao) for the rubidium dimer. The decay of the ratio at large 
r is fastest for large n. 



be introduced in pairs. There is no point in including 
the Vn potential without also including the V12 poten- 
tial. Indeed, inclusion of just one member of the (Vn,Vi2) 
would most likely degrade rather than improve the ac- 
curacy of the dispersion potential. Given that the Cq 
parameter has been calculated to a precision of better 
than 1% for most alkali systems [2lL l24l |25| . usage of a 
dispersion interaction involving the (Vn, V12) potentials 
may be warranted. 

There have been a couple of experimental investiga- 
tions of alkali-dimer potentials that have included dis- 
persion forces with n > 10. The value of C\\ for the 
Rb dimer has been given by van Kempen et al . How- 
ever, this value of C\\ is most likely model dependent 



for reasons discussed earlier. Consideration of the V\i(r) 
potential has also occurred in analyses of the spectrum 
of the Cs dimer [sSE3- In this case, the Vu(r) poten- 
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FIG. 6: The ratio of |K(r)|/|Wi 6 (r)| for n = 6, 8 and 10 as 
a function of r (units of ao) for the sodium dimer. The other 
curves show the ratio of |V2ra-i(r) + V2n(r)|/|Wi6(r)| for 2n = 
12, 14, 16. 
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FIG. 7: The ratio of \V„(r)\/\W 16 (r)\ for n = 6, 8 and 10 
as a function of r (units of ao) for the rubidium dimer. The 
other curves show the ratio of \V2n~i{r) + V2 n (r)\/\Wi6(r)\ 
for 2n = 12, 14, 16. 



tial was omitted so it is doubtful whether the inclusion 
of Vi2(r) was justified in this case. 
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